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\q 1 Abstract 
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Utilizing an appropriate ansatz to the wave function, we reproduce the exact 

in 

fT) ■ bound-state solutions of the radial Schrodinger equation to various exactly 

solvable sextic anharmonic oscillator and confining perturbed Coulomb mod- 
. els in D-dimensions. We show that the perturbed Coulomb problem with 



a: 



eigenvalue E can be transformed to a sextic anharmonic oscillator problem 



with eigenvalue E. We also check the explicit relevance of these two related 

& . 

problems in higher-space dimensions. It is shown that exact solutions of these 
potentials exist when their coupling parameters with k = D + 21 appearing 
in the wave equation satisfy certain constraints. 
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I. INTRODUCTION 



The solution of the fundamental dynamical equations is an interesting phenomenon be- 
cause of its importance in quantum-field theory, molecular physics, solid-state and statistical 
physics. To obtain the exact £-state solutions of the Schrodinger equation (SE) to various 
quantum mechanical problems are possible only for few potentials and hence approximation 
methods are used to obtain their solutions. According to the Schrodinger formulation of 
quantum mechanics, a total wave function provides implicitly all relevant information about 
the behaviour of a physical system. Hence if it is exactly solvable for a given potential, 
the wave function can describe such a system completely. Until now, many efforts have 
been made to solve the stationary SE with a sextic anharmonic oscillator and perturbed 
Coulomb potentials in one to three dimensions through the Hill determinant matrix method 
[1-4]. The study of the SE with these potentials provides us with insight into the physical 
problem under consideration. Further, the study of SE with some of these potentials in the 
arbitrary dimensions D is presented [3]. 

The purpose of this paper is to carry out the analytical solutions of the D-dimensional 
radial SE with exactly-solvable Coulomb plus linear plus harmonic (CLH) V\{r) = —a/r + 
br + cr 2 and sextic anharmonic oscillator (AHO) V 2 (r) = fir 2 + Ar 4 + r/r 6 potentials through 
an appropriate ansatz to the wave function. In cases of exactly solvable models, the wave 
function can be expressed as a finite power series polynomial multiplied by an appropriate 
reference function (usually, the asymptotic form) to reproduce the exact solutions. The 
analytical solution of the Schrodinger equation for the energy levels with a class of confining 
potentials of type V\ have been studied by Datta and Mukherjee [5] by using Kato-Rellich 
perturbation theory for linear operations. It is well known that this confinement potential 
has been used for calculation of qq bound-state masses [6]. Killingbeck [7] has calculated 
the energy eigenvalues of the confinement potential by using hypervirial relations. Exact 
solutions [8] of potentials of type V\ and V2 are obtained by a number of authors in three- 
dimensional space when the coupling parameters satisfy certain relations. Chaudhuri and 



2 



Mondal [3] studied the D-dimensional sextic AHO and CLH problems within the frame- 
work of supersymmetric quantum mechanics (SUSYQM) and Hill determinant method [4,9- 
11]. Chaudhuri and Mondal have shown that SUSYQM yields exact solutions for a single 
state only for the quasi-exactly-solvable potentials of type V\ and V 2 in D-dimensions with 
some constraints on the coupling constants [3]. They have also obtained numerical results 
throughout the Hill determinant method. The ideas of supersymmetric quantum mechanics 
have been used for the study of atomic systems [11], the evaluation of the eigenvalues of 
a bistable potential [12], the improvement of the large-N expansion [13], the analysis of all 
known shape invariant potentials [10,14], and the development of a more accurate WKB 
approximation [14]. Tymczak et al [4] devised a highly accurate quantization procedures 
for the inner product representation both in configuration and momentum space for vari- 
ous AHO potentials in one and two dimensions. Additionally, Dobrovolska and Tutik [15] 
studied the bound-state problem within the framework of the SE through the logarithmic 
perturbation theory. Recently, they also extended the formalism to the bound-state problem 
for spherical oscillator of type V 2 with its subsequent application to the doubly anharmonic 
oscillator [16]. Furthermore, a simple formalism [17] based on a suitable choice of the wave 
function ansatz has been proposed for reproducing exact bound-state energy eigenvalues and 
eigenfunctions for exactly solvable model within the framework of the SE. Very recently, this 
simple approach has also been applied [18] with remarkable success, to various molecular 
quantum mechanical problems in D-dimensions [19,20]. 

The object of this paper is to extend the above simple approach to reproduce bound- 
state exact solutions for potentials of type V\ and V 2 in D-dimensions with some constraints 
on the coupling constants. For a certain choice of parameters the method provides exactly- 
solvable potentials for the D-dimensional sextic AHO and CLH problems. We then compare 
our results with the exact ones obtained from the transformation of LHO into sextic AHO 
problem. 

This paper is organized as follows. In Section II, we solve analytically the D-dimensional 
radial Schrodinger equation for the sextic AHO and CLH problems by a suitable choice of a 
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wave function ansatz to each exactly-solvable problem. On the other hand, the exact energy 
eigenvalues of sextic AHO are obtained by transforming radial wave SE in D-dimensions with 
angular momentum i to another problem in (2D — 4)-dimensions with angular momentum 
2£ + 1. The results and conclusion will be given in Section III. 

II. THE D-DIMENSIONAL RADIAL SCHRODINGER EQUATION 

In the D-dimensional Hilbert space, the reduced radial wave Schrodinger equation (with 
h — m—1 units) for a spherically symmetric potential V(r) takes the form [21] 



rfl + (^l)rf_«i±B-2) +2(£ _ y(r)) 



dr 



r dr 



V>(r) = 0, 



(1) 



where the interaction potential is chosen to be of type V\ or V 2 and E stands for its eigen- 
values. Further, equation (1) can be simply transformed to the form [21] 

f d 2 [(k — l)(k — 3)] 



dr 2 



4^2 



+ 2 (E - V(r)) }R(r) = 0, 



(2) 



where R(r), the reduced radial wave function, is defined by 



i?( r ) = r P-i)/2^( r ) > 



(3) 



and 



k = D + 2£, 



(4) 



which is a parameter depends on a linear combination of the spatial dimensions D and the 
angular momentum quantum number I [21]. 

We substitute r = 7p 2 /2 and ip = x(p)/p to transform Eq.(l) into another Schrodinger- 
type equation in (D' = 2D — 4)-dimensional space with angular momentum L — 2£ + 1, 



p dp 



where 
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E = ^p 2 E, V(p) = 7 V\/( 7 p 2 /2), 7 = V(-^) 1/2 . (6) 

It is seen that through this transformation, the D-dimensional radial wave Schrddinger 
equation (1) with angular momentum £ can be transformed to a (D' = 2D — 4)-dimensional 
problem with new angular momentum L — 2£ + 1. In particular, under this transformation, 
CLH problem of type V± with eigenvalue E can be transformed to a sextic AHO problem of 
type V2 with eigenvalue E as 

t/(p)=/i P 2 + Ap 4 + ^ 6 , (7) 
with coupling constants given by 



and eigenvalue 



A. Confining perturbed Coulomb problem 

We attempt to solve the wave equation (2) of reduced radial wave R(r) in the IP- 
dimensions for a spherically symmetric potential of confining CLH form [3,22]: 

V(r) = -- + br + cr 2 , c> 0. (10) 
r 

In particular, for c = and b > 0, such a potential reduces to the well known quarkonium 
Cornell potential (cf. Refs. [21] and references therein). Apart from its releavance in heavy 
quarkonium spectroscopy (cf [21] and references therein), this class of potentials with c = 
has important applications in atomic physics. For exactly solvable problems such as CLH, 
the representation of the radial portion of wave function ansatz, containing an appropriate 
reference function (usually, asymptotic form), is 

R(r) = exp [p (a, [3, r)] £ a n r 2n+ ^'\ (11) 

n=0 
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provided that the power of the reference function has the following selection: 

p (a, [3, r) = -ar 2 + (3r, 



(12) 



should fall faster than the asymptotic form of the wave function. Upon substituting Eq. 
(11) into Eq. (2) and equating the coefficients of r 2n+ ( fe ~ 1 )/ 2 to zero, we readily arrive at the 
following relation 



A n a n + B n+ ia n+ i + C n+ 2 a n+2 — 



(13) 



where 



A n = 2E + f3 2 + a (An + k) , 
B n = 2a + (3 (An + k - 1) , 

C n = An 2 + 2n(k-2), 

and the value of the parameters for p (a, f3, r) can be evaluated as 

b 

a = ±^2~c, /3 = ± 



(14) 



(15) 



(16) 



2c 



(17) 



To obtain a well-behaved solution at the origin and infinity, it is more convenient to take 
a = —\[2c and (3 = — -j= which ensure that wave function ansatz representation in (11), 
be finite for all r and convergent at oo. Further, for a given p, if a p ^ 0, but a p+ i = a p+2 = 
a p+3 = ■ ■ ■ = 0, we then obtain A p = from Eq. (14), i.e., 

(18) 



= -- + x /-(4p + 2£ + D), p = 0,l,2, 



Carrying through a parallel analysis to Refs [17,18], A n , B n and C n must satisfy the deter- 
minant relation for a nontrivial solution 

Bo 

A B 1 C 2 



Det 



A p _ x B p 



0. 



(19) 
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To utilize this method showing the simplicity of this approach, we pursue determinant 
analysis to present the exact solution for p = 0, 1 as follows. 

Case (1): when p = 0, we get from Eq. (18), the exact bound-state solution of the CLH 
problem in D-dimensions. So that the eigenvalue (ground state) is given by [17,18] 



E° = — + f-{2l + D) . 



(20) 



Further, it is shown from Eq. (19) that B = 0, which leads to the following constraint on 
the coupling parameters as 



b(k - 1) = labile. 



(21) 



which consequently, from Eqs. (20) and (21), the perturbed Coulomb potential admits an 
exact ground eigen energy: 



b(k-\f b(k-l) 



4a 2 



2a 



2a 



(22) 



(k-1) 2 

which is consistent with Refs [3,22]. Further, the corresponding wave function (unormal- 
ized): 

2ar b(k — l)r 2 



4 fc) (r) = a /exp L (fc _ ^ ^ 
Case (2): When p — 1, the exact energy spectrum becomes 



(23) 



E« = - 
1 2 



6(A; + 1) 2 b(k + l) 



Aa 2 



(24) 



(25) 



2a 2a (k + 1) 2 _ 

and the corresponding wave function (unnormalized) can be readily found as: 

(k)/ t 2ar b(k + l)r 2 
y>i{r) = (a + a 1 r)rexp Ta ' 

where ao and a\ are normalization constants. The relation between them can be determined 
through the relation B a + C\d\ = to be 

b la 



ai = 



V2~c (k-1) 



a . 



(26) 
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Furthermore, we have the following recurrence relation from Eq. (19), that is, B B 1 = A Q C\, 

which consequently provides the following constraint on the coupling constants of potential: 

b 2 



4a 2 - "-(k - l)(k + I) 



^Kk = ^k(k-l). 



(27) 



Following this approach, we can further generate a class of exact solutions through setting 
p — 0, 1, 2, • • • , etc. Generally speaking, if a p ^ 0, a p+ i = a p+ 2 = • • • = 0, from which we can 
obtain the energy spectra (cf. determinant (19)). For the generalization, one needs to use 
the shape invariance property and the relation between supersymmetric partners [22,23] to 
find the general solution of energy eigenvalues as 



El 



b(2n + k-lf + b(2n + k-l) 



4o 2 



2a 



2a 



with the corresponding wave functions 



(p), 



(ao + air + • • • + a p r p ) r exp 



2ar 



(2n + k- l) 2 

b(2n + k-l)r 2 



4a 



(28) 



(29) 



(2n + k- 1) 

where a^ (i — 0, 1, 2, • • • ,p) are normalization constants. 

Finally, considering the following exactly solvable potentials. In particular: (i) Harmonic 
oscillator: when a = b = and c = \uo 2 are inserted to Eq. (10), giving a = —uo and (3 = 0. 
Thus, we can readily obtain the energy eigenvalues through using Eqs. (14)-(16) as [16,18,24] 



E n£ = - (An + D + 21) , n,^ = 0,l,2, 



(30) 



(31) 



and the corresponding radial wave function becomes 

^W(r) = (a + air 2 H h a„r 2n ) / exp 

where ai with i — 0, 1, 2, • • • , n are normalization constants. 

(ii) Coulomb problem: when b = c = 0, and a — Z, implies that a = and f3 = 



• (2n+D+2i-i) ' f rom which we readily obtain the exact eigenvalues as [25] 

2Z 2 



nl ~ (2n + D + 2£ - 1) 
and radial wave function becomes 



2 - 



n,£ = 0,1,2, 



^W(r) = (a + a^ 2 + ■■■ + a n r 2n ) / exp 



2Z 



(2n + D + 2£ - 1) 
where a^ with i — 0, 1, 2, • • • , n are normalization constants 



(32) 



(33) 
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B. The sextic anharmonic oscillator problem 



This version of sextic AHO potential 

V(r) = fir 2 + Ar 4 + r/r 6 , 77 > 0, (34) 

has been studied in the D dimensions through Hill determinant method [3]. We want to 
solve the radial SE, Eq. (2), with Eq. (34) by selecting the following representation of 
ansatz to the radial portion of wave function 

R(r) = exp [p (a, [3, r)] ]T a n r 2n+ ^'\ (35) 

n=0 

provided that the power of the reference function has the following selection: 

p (a, [3, r) = -ar 2 + -flr 4 . (36) 
2 4 

Implementing the present method on the representation (35), and taking the coefficients of 
r 2n+(fc+i)/2 ^ zer0; we obtain the relation 

A n a n + B n+ ia n+ i + C n+ 2d n+ 2 = 0, (37) 

where 

A n = a 2 + (3 {An + k + 2) - 2/x, (38) 

B n = 2E + a (An + k) , (39) 

C n = 2n(2n + k-2), (40) 
and the value of the parameters for p (a, j3, r) can be evaluated as 

P = ±j2~n, a = ±^=. (41) 



To obtain a well-behaved solution at the origin and infinity, we must set (3 = —^/2rj and 
a = —^Fj wn i cn ensures that wave function ansatz, Eq. (35), be finite for all r. Further, 
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for a given p, if a p ^ 0, but a p+ \ = a p+2 = a p+3 = ■ ■ ■ = 0, we then obtain A p = from 
Eq. (38), which leads to the following constraint on the coupling parameters of sextic AHO 
problem as 



2/i + J^q{An + k + 2) - ^- = 0. 
v 2?] 



(42) 



Case (1): when p — 0, it is shown from Eq. (19) that _£>o = 0, which leads to the following 
energy eigenvalue (ground state): 

Xk 



rpk 

" 2^' 

and the corresponding wave function (unormalized) : 

Ar 2 + r\r A 



(43) 



ijjo{r) = a r exp 



v 7 ^ 



(44) 



Case (2): when p — 1, it is shown from Eq. (19) that -£>o-E>i = A Ci, which leads to the 
constraints on the coupling parameters of the potential as 

A 



El 



'A 2 



A: 



2^ (fc+2)+ V5; ( ' c+2) -2 



\[2r)(k + 2) + 2/i 



(45) 



In particular, an important version of the above sextic AHO is the harmonic oscillator 
problem 



V(r) = fir 2 , n > 0. 



Selecting the following representation of the wave function ansatz:: 



(46) 



(47) 



with 



n=0 



1 



(48) 



and itirating the previous steps, one gets: 



A n = a 2 - 2/i, 



(49) 
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B n = 2E + a (An + k) , 



(50) 



C n = 2n(2n + k-2), 



(51) 



which leads to choosing a = —yfljl. 

Case (1): when p — 0, it is shown from Eq. (19) that B = 0, which leads to the following 
energy eigenvalue (ground state): 

^/2j2k 



En 



and the corresponding wave function (unormalized) : 



i/>o k \r) = a r*exp 



(52) 



2\i 2 

r 



(53) 



Case (2): when p — 1, it is shown from Eq. (19) that B B 1 = A Ci, which leads to the 
restriction on k and the parameters of the potential as 



El 



(k + 2) + J2fi. 



(54) 



Generally speaking, the energy eigenvalues are 

v/2/i 



Et 



(An + k), 



(55) 



and the corresponding wave function can be read 

■0n( r ) — (°o + ai r H 1- <V 2 ™) ^ exp 



2/* 2 

r 



(56) 



where (i = 0, 1, 2, • • • , n) are normalization constants. 



III. CONCLUDING REMARKS 

We applied the wavefunction ansatz method to the confining CLH and sextic AHO 
interactions. Table 1 shows the calculated energies of the CLH type potential together with 
those obtained by the SUSYQM and Hill determinant method for high values of parameters 

11 



in three- and four- dimensions. We know from (8) and (9) that the CLH problem with exactly 
solvable potentials Vi^'^'^ 3 \r) in four- dimensions can be transformed to the sexrtic AHO 
problem (7) in four- dimensions with the exact eigenvalues given through (9). 

We compute the energy eigenvalues of the following conjugate AHO potentials: 



\3/2 
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V/ 4, (r) = r 2 + [l/2(7.625) 3 / 2 J r 4 + [l/4(32)(7.625) 2 ] r 6 , Vf j (r) = r 2 + [1/2(7.375) 
+ [l/4(32)(7.375) 2 ]r 6 and vl %) {r) = r 2 + [l/2(7.125) 3 / 2 ] r 4 + [1/4(32)(7.125) 2 ] r 6 in two- 
dimensions by the present simple method and compute our results in Table 2 with the 
exact values given by (9). A class of AHO may be constructed from (43) that ad- 
mits exact solutions. We also compute the energy eigenvalues of the following conju- 
gate AHO potentials: t// 4) (r) = r 2 + [l/2(7.5) 3 / 2 ] r 4 + [l/4(32)(7.5) 2 ] r 6 , v} 5) (r) = r 2 + 
l/2(7.25) 3 / 2 ] r 4 + [1/4(32) (7.25) 2 ] r 6 and t// 6) (r) = r 2 + [l/2(7.0) 3 / 2 ] r 4 + [1/4(32) (7.0) 2 ] r 6 
in four-dimensions by the method and compute our results in Table 3 with the exact values 
given by (9). These eigenvalues are checked by other methods. 

This method yields exact solutions for a single state only for a potential of type V 2 and 
many states of type V 2 with some constraints on the coupling parameters. Our method 
is applicable to many general CLH or AHO and produces excellent results for the low- 
lying states. It gives the exact solutions of the Couloimb and the harmonic oscillator in 
.D-dimensions. The CLH and sextic AHO in D-dimensions are related through Eq.(9) and 
are verified in Table 2 by this method. 

A class of conjugate AHO having exact eigenvalues may be constructed from the trans- 
formation of CLH potential. It is found that the eigenvalues of central potential V(r) are 
identical for D = 2,1 = 2, D = 4, £ = 1 and D = 6, £ = states. This is because k = D + 2£ 
remains unaltered for these states. 
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TABLES 

TABLE I. The lowest energy eigenvalues of the CLH potential for I = 0,1,2 in three- and 
four-dimensions . 



a a 


b 


c 


i 


D 


Present method 


SUSYQM [3] 



12 
6 

10 
14 



32 

J_ 

32 

J_ 

32 

J_ 

32 

J_ 

32 




1 

2 

1 
2 



3 
3 
3 
4 
4 
4 



-7.625 
-7.375 
-7.125 
-7.500 
-7.250 
-7.000 



-7.625 
-7.375 
-7.125 
-7.500 
-7.250 
-7.000 



a The parameter values here are taken from [3]. 

TABLE II. The eigenvalues of the conjugate sextic anharmonic oscillators in two-dimen- 
sions. are compared with the exact values. 



Conjugate sextic AHO a 



Present work 



Exact Value, Eq.(9) 



V} 4 \r) 



2.8971438733606 
5.8916775545493 
8.9912237911843 



2.8971438733606 
5.8916775545493 
8.9912237911846 



i The parameter values in constructing the sextic AHO are taken from Table 1. 



TABLE III. The eigenvalues of the conjugate sextic anharmonic oscillators in four-dimen- 
sions. are compared with the other works and the exact values. 



Conjugate sextic AHO a 



Present work 



Hill Determinant [3] Exact Value, Eq.(9) 



V} 4) (r) 



1 4.3817804600412 
3 7.427813527082 
5 10.583005244257 



4.381780461 
7.427813527 
10.583005244 



4.381780459 
7.427813526 
10.583005240 



A The parameter values in constructing the sextic AHO are taken from Table 1. 
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